Formulation of Time-Resolved Counting Statistics Based on a 
Positive-Operator- Valued Measure 
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We propose a derivation of the full counting statistics of electronic current based on a positive- 
operator-valued measure. Our approach justifies the Levitov-Lesovik formula in the long-time limit, 
but can be generalized to the detection of finite-frequency noise correlations. The combined action 
of the projection postulate and the quantum formula for current noise at high frequencies imply an 
additional white noise. Estimates for this additional noise are in accordance with known experi- 
ments. We propose an experimental test of our conjecture by a simultaneous measurement of high- 
and low-frequency noise. 



The core of quantum measurement theory is the pro- 
jection postulate p|. It provides a consistent description 
of a sequence of measurements. Quantities represented 
by non-commuting operators cannot be measured simul- 
taneously. The corresponding projection operators have 
to be time-ordered. For continuous variables the projec- 
tion postulate should be replaced by a positive-operator- 
valued measure (POVM) [1]. The idea of the POVM is 
that one does not measure the exact value for a given op- 
erator but a finite accuracy is taken into account due to 
some interaction with the detector and its internal dy- 
namics. However, due to Naimark's theorem Q, every 
POVM can be realized by a set of orthogonal projections 
in an extended Hilbert space. The resulting POVM will 
depend on the detection scheme. 

The statistical behavior of current flow in a quantum 
point contact can be found by measurement of the corre- 
lation functions. The long-time cumulants of the trans- 
ferred charge can be derived from the Levitov-Lesovik 
formula Q , which led to the foundation of the electronic 
version of full counting statistics (PCS) 0]- It has 

been confirmed experimentally for noise [8| and third 
cumulant I, [13, El. o n the other hand, the current- 
current correlation function (noise spectral density) is 
given by the quantum noise l^i l3] j which coincides with 
the PCS result at low frequencies. The high-frequency 
quantum noise can be obtained by a generalization of 
PCS to finite frequencies with additional predictions for 
higher cumulants at arbitrary frequency lj,ll5| . Also the 
semiclassical predictions of the third cumulant are con- 
sistent with purely quantum results in some limits [l^ . 
The behavior of quantum noise has been confirmed ex- 
perimentally also for high frequencies 17, 18, 3|. Prom 
the fundamental point of view, the low-frequency results 
can be justified by a proper use of the projection postu- 
late but there is no unique derivation for high frequencies 
[20| . A similar problem occurs for a chain of spin-resolved 
detectors, for which the results depend on the detector 
properties ^ . While it is reasonable to expect an influ- 
ence of the detector on the outcome 0,[21, 22 1, it should 



be possible to separate it from the bare signals of the 
sample. 

In this Letter we address the question, if a general defi- 
nition of PCS for finite frequencies is possible - maintain- 
ing the probabilistic interpretation. We will demonstrate 
that the standard definition of PCS, when generalized to 
finite frequency, can lead to negative probabilities. To 
cure this deficiency, we show that taking into account a 
minimal model of a detector, a POVM of PCS can be 
introduced, which leads to positive definite probabilities. 

The definition of the generating functional for a prob- 
ability distribution p of a given time trace of the current 
through a quantum point contact, is 



ix{t)I{t)dt/e 



= JdI g[I]ef 



ix{t)I{t)dt/e 



(1) 



On the other hand, one can first define PCS generating 
function 0,1, 

„ (2) 

Here p denotes initial state density matrix, I(t) is the 
Heisenberg current operator, (j> refers to classical phase 
bias, and T(T) denotes (anti-)time ordering. A detailed 
definition of I{t) will be given later. 

Taking S[x] = S[x, 0] we obtain g by inverse Pourier 
transform of ([1]). However, this gives positive probabili- 
ties only in the zero- frequency limit. Por time-dependent 
quantities, we can construct the following counterexam- 
ple for a single mode point contact at zero temperature in 
the tunneling limit (transmission T <C 1). Let us define 



X = 



ta 



dtdt'I{t)Iit')[e-^'-'"^"^'" -2e-^'-''^"/^'"]. (3) 



Then, following [ij we find ((^X)^)^ = -Ttoe^/SsTr^/^ 
for 5X = X — {X)g and to ^ s. This obviously contra- 
dicts the interpretation of g as probability. 

To overcome this fundamental problem, we now con- 
struct a positive definite probability of time-dependent 



2 



FCS based on a POVM. Instead of the projection oper- 
ator we define the more general Kraus operator (23j 



k[I] = / DipTe^'^''^ 



tvit)lHt)~Iit)]/e-^\t)/r 



(4) 



Causality is preserved since the detector affects the mea- 
surement only in later times. The time scale r describes 
internal fluctuations of the detector and depends on its 
temperature in general. For r — > cx) the measurement is 
accurate, but the detector noise strongly affects the sys- 
tem by full projection. A shorter r reduces the influence 
of detector but induces a larger measurement error. The 
integration measure contains also a normalization factor 
to be determined later. The positive definite probability 
of a given I{t) is defined as 



p[I]^TrpK^[I]K[I], 



(5) 



for the given initial density matrix p. We note that our 
choice of the Kraus operator represents generically the 
influence of a detector, parametrized by a single param- 
eter T. What concrete models of detectors lead to our 
definition of the Kraus operator is an interesting ques- 
tion, which we will not address here. 

We now substitute in Eq. ([5]) (p 4> zk x/'^ in K 
and , respectively. The generating functional S[x] — 
ln(exp(i J dtx{t) I (t) / e)) p necessary for the calculation 
cumulants takes the form 



=ln J D<j>> 



,S[x.0I-/dt[20=(t)+x"(t)/2]/T 



(6) 



where S[x, 0] is defined by The measure Dcj) is scaled 
to keep S[x = 0] = 0. The measuring device affects 
the generating function by the additional exponent in 
([6]). In 121]] Di Lorenzo and Nazarov used the expression 
t"^ ^"^ + cjp / {Acj})"^ instead of (j)'^, with Acj) as an additional 
parameter, and considered low-frequency measurements. 
In contrast, we rather assume a continuous weak mea- 
surement of the system to obtain finite frequency corre- 
lations. 

To further model our measuring device, we note that 
in general a current measurement has also a spatial sen- 
sitivity, i.e. the point of the measurement is not exact. 
In experiments, it can be usually related to the finite ca- 
pacitance of the sample. Therefore, we assume a generic 
form of the current operator in a quasi-one-dimensional 
lead as 



m 



i{x, t)e~ *"2A?2 dx/{V2T:Ax) . (7) 



The setup is shown in Fig. [T^. The real dispersion may 
be non-Gaussian. However, we stress that our model is 
general enough to capture the essential physics, but still 
allows some analytical progress. 

We will assume non-interacting electrons and energy- 
and spin-independent transmission through the Ad mode 




FIG. 1: (a) Sensitivity of current measurement. The Gaussian 
distribution refers to the dispersion of current measurement, 
(b) The function q{z). 



junction. We count all modes, although most of them are 
just reflected and denote the Fermi velocity, the trans- 
mission and the reflection coefficients for mode n by d„, 
T„ and i?„ — 1 — T„, respectively. For convenience, we 
introduce t„ = |a;o|/i^n and t„ = Ax/vn- The times r„ 
are related to RC times of the circuit, which limits the 
observable frequencies to a; < t~^. Furthermore, we as- 
sume that <C T~^, which means that the detector 
sensitivity function is entirely located on one side of the 
junction. 

To model the electron transport we apply the stan- 
dard scattering picture around the Fermi level [l^ . 
The Hamiltonian can be approximated by H = 
J2nl dxTin{x), where 

Ufl = ihVn[lp\^,--^{x)dx4'Ln{x) - ■lpl^f^{x)dx-4'B.n(x)] (8) 
+ qn6{x)['4!lf^{x)lpRf,{x) + ^Rf,ix)4>Lnix)] 

-eV0ix)[i,i^{x)i^Ln{x) + i^l^ix)i^m{x)]. 

The scattering states obey standard fermionic anticom- 
mutation relations {■4'\n{x),ipBm{y)} = SABSnm6{x - y) 
and {i'An{x),^Bm{y)} — 0. Here A = L,R denote left 
and right going state, n = (n, cr) denotes mode number 
n and spin orientation a. The transmission coefficient 
is given by r„ — cosh^^((7„/?iu„). The current operator 
is defined as /(x) = J2n'^'"nipLfii^)'^Ln{x) - L ^ R. 
The initial density matrix for a thermal state is p = 
g-H/ksT ^rY^^-H/kBT time evolution is governed 

by the Heisenberg operator I{x, t) = e^^*'/^I{x)e~'^^*'/^ . 

For our model, the mean current is independent of 
the detector, {I{t))p = GV, where the conductance 
G = En^nGg and Gq = e'^/irh. We define the 
noise spectral density as a second cumulant e^Piuj) — 
J dte^^^{5I{Q)5I{t))p, where 51{t) = I{t) - {I{t))p. It is 
calculated from the functional derivative 



(«(0)«(0)p = 



5'S[x] 



Sx{t)Sx{0) 



(9) 



where S[x] is defined by Eq. In our construction, the 
noise is a classical quantity and, hence, symmetric with 
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respect to w. We consider frequencies |wr„| <C 1, since 
we do not include capacitive effects and obtain 



P{u;) = Pos + Psiuj) + Fo(^) + Pr{cj) + PA{io)- 



(10) 



Let us discuss the behavior of all terms of this expres- 
sion. The first one, Poff = 1/t is a white offset noise, 
independent of temperature and voltage bias. Defining 
w{u}) = u}cth{huj /2kBT) and w±{llj) = w{llj ± eV/h) the 
second term 

Ps(^) = 5]^{2r„u;(a;)+i?„[w+(c^)+zi;_(c^)]} (11) 



is just the symmetrized quantum noise 
J dtcos{ujt)TTpSi{0)6i{t)/e'^ [13]. However, for 

energy independent transmission, as we assume 
here, the asymmetric noise contains only the ad- 
ditional term X^n^"^/""' which is independent 
of temperature and voltage. The next term is 
Po{lj) = ^„ 2i?„ sin^(a;t„)w(aj)/7r. This is a con- 
tribution to the quantum noise due to the finite flight 
time to the detector, as it depends on i„. Note that it is 
independent of voltage and sensitivity. The problem of 
flight time has been already discussed in context of third 
cumulant l^, 24 1 , but there is no experimental evidence 
of its influence on the noise. The detection noise, 



PAuj) 



T 

4^ 



E 



Uj{l + RnC 



2iLjtn 



(12) 



combines the effects of the measurement sensitivity t and 
flight times tn but is independent of voltage and temper- 
ature. Finally, 

^fni^- a)RnTn [w+ {a)+W- {a)] 

(13) 

is an additional mixed noise. Here the sensitivity ampli- 
tude is given by 



fn{a) 



dt 



exp 



i-e-*V4-r 

80rT„/T 



(14) 



It is independent of the flight times, but all other pa- 
rameters enter in a rather complicated way. Eq. is 
the only term depending on r„ in the limit \ujTn\ <C 1 
However, voltage and temperature are arbitrary. For 
T <C Tn the amplitude reduces to /n(a) = T(e~" — 
,5(a)0F/T„)/4. 

We assume that most modes are closed, ^„ T„ <C M, 
which is true in many experimental setups (e. g., tunnel 
barrier, diffusive wire, quantum point contact). We will 
consider several interesting limits: short and long wire 
(flight time) , low and high frequencies and zero tempera- 
ture. They correspond to the most common experimen- 
tal setups. For short flight times, |w<„| <C 1, we have 



Po = 2a;^t(;(w) X]„i^- For long flight times, \ujtn\ > 1, 
we get Pq = Mw{ijj) /tt because random flight times imply 
sin2(wi„) 1/2. In both cases P^ = (E„ "^n ^)^'^/47r^ 
is independent of voltage and yields additional offset 
noise. For low frequency and a slow detector, h\uj\ 
{h/r, h/Tn) ^ eV, UbT the mixed noise is negligible since 
Pa < Ps(0). However, for ksT < h\uj\ < fi/^n < ^/t, 
we have Pa = (t/Stt^) P„Png(r„|eF|/S)/r2, ^j^^^g 



2z e * dt. The decay of q{z) is shown 



in FigdlD. It vanishes for eV 3> h/Tn, which means that 
the size of the wave packet becomes smaller than the 
spatial sensitivity of the detector. 

From the above results we conclude, that the POVM 
reproduces the standard quantum result in the case \uj\ <^ 
(l/r„,l/<„) < 1/t < eV/h,kBT/h. It can be shown, 
that corrections to higher zero-frequency cumulants are 
negligible in this case, because Pa is small. Hence, the 
Levitov-Lesovik formula and FCS [J,|^,0,[3 are justified, 
as expected. 

The situation becomes more interesting, if we look at 
the high frequency quantum noise. We can make P^- and 
Pa negligible by choosing a very small r, which corre- 
sponds to a weak detection. For small tn, also Pq gives 
only a small contribution. Moreover, higher cumulants 
then have also negligible corrections to predictions of gen- 
eralized FCS as small r corresponds to (/i ^ in eq. ([2]). 
What remains is the large white Gaussian offset noise Pog 
- the price we have to pay for small t. Conversely, lower- 
ing Poff will increase Pa, which additionally depends on 
voltage. They become of the same order at t ~ t„ and 
Pa is growing as yr/Tn for r ^ t„. Hence, one cannot 
get rid of the additional noise by increasing t, since it 
increases the backaction noise. 

The offset noise for a few mode quantum point contact 
in most high frequency experiments 17|, |18| was usually 
subtracted. However, the results of the recent experi- 
ments [3 [lo'l show relatively high absolute noise temper- 
ature T* = CqPh/^GkB, setting an upper bound to the 
offset-noise. We find T* = Tos + Ta with Toff = 24K and 
Ta — 1.5mK for a single-mode quantum point contact 
with transmission 1/2, r = Ips and ti = lOps. Although 
Toff is larger, it is constant whereas Ta depends on volt- 
age bias and drops to zero according to the function q{z). 
For zero temperature this yields a characteristic voltage 
of 130/iV at z = 1. Diffusive and tunnel barriers have 
usually a higher conductance due to the large number 
of modes and hence the offset temperature can be much 
lower. On the other hand, in a Josephson junction or 
a quantum dot as a detector, the measured quantity is 
more qualitatively than quantitatively related to the fre- 
quency dependent quantum noise [25 1. 

Moreover, as the offset noise is white, it must be visi- 
ble also at low frequencies. It would be therefore of great 
interest to test high and low frequency noise simultane- 
ously. The high frequency detector measures only the 
difference P{eV) — P(0) whereas at low frequencies the 
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absolute noise is measured. The fluctuation-dissipation 
theorem [2^ will of course be maintained after subtrac- 
tion of the offset noise and then taking the limit t ^ 0. 

We generalize our approach to a multi-terminal mea- 
surement. We define the Kraus operator ^ as 

K[I] = j D^re-^''*^'*^*'^-***)!^"'*^*)"-^^^*)!/''"'^^^*'/^-*^ 

(15) 

where A labels the terminals. A simple example is 
the noise between the left and right sides of a junc- 
tion, namely e'^PABi^^) ^ J dt e^'^\5lA{Q)5lB{t)) where 
A, B = L, R for left and right terminals (or equivalently 
xq < and xo > 0), respectively. The cross correlation 
Plr{^) is finite only for |u;t„| <C 1, since otherwise one 
averages over different t„. In this limit, we have 

PlL = PloS + Ps + Prr. + PaL + PaR, 
PrR = ProS + Ps + Prn + PaL + PaR, 

Plr = Prl = -Ps + P.- Pal - Par, (16) 

where P^{uj) = uj^tl + tr)M ^„ Tjiir^ Here P^ffA = 
T^^, Pr^ and Paa are given by Eqs. (fT2|) and (fT3|) with 
T and T„ replaced by ta and t„^. Pa is here replaced 
by the sum of contributions of both detectors. The cross 
noise does not contain the offset noise, so it can help 
to estimate the measurement timescales, in particular 
Pal + Par at low frequency. 

Finally, we propose the following test of our definition 
of a quantum probability. Consider two detectors, similar 
to [10[, measure Xl and Xj^^ defined by ([3]) for J = II 

cl 

and In, respectively. The classical expectation Xl = Xji 
based on charge conservation in the low frequency limit, 
leads to {SX^SXiijci > 0. The quantum measurement 
using the probability density p[I], results for r ^ s in 
{SXlSXj^} p = {{6X)'^)g. This is, however, negative as 
we have shown in beginning. 

In conclusion, we have constructed a positive proba- 
bility measure, based on POVM, that justifies the use of 
Levitov-Lesovik formula and FCS in long time (low fre- 
quency) limit. Our approach cures certain deficits in the 
standard definition of FCS, which lead to negative prob- 
abilities. Introducing a generic one-parameter model of 
the influence of the detector, we predict an intrinsic ad- 
ditional white offset noise. Such a noise is in agreement 
with recent experiments and a further verification by si- 
multaneous low- and high-frequency noise measurements 
would be desirable. 
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